TRAVELING WAVE DISPERSAL IN PARTIALLY SEDENTARY 
AGE-STRUCTURED POPULATIONS 



THUC MANH LE, FRITHJOF LUTSCHER, AND NGUYEN VAN MINH 



Abstract. In this paper we present a thorough study on the existence of 
traveling waves in a mathematical model of dispersal in a partially sedentary 
age-structured population. This type of model was first proposed by Veit and 
Lewis in [Am. Nat, 148 (1996), 255-274]. We choose the fecundity function to 
be the Beverton-Holt type function. We extend the theory of traveling waves 
in the population genetics model of Weinberger in [SIAM J. Math. Anal, 
13 (1982), 353-396] to the case when migration depends on age groups and a 
fraction of the population does not migrate. 



1. Introduction 

In [24] a mathematical model for dispersal in a partially, sedentary age-structured 
population was developed to simulate the spatial spread of the house finch (car- 
podacus mexicanus). The house finch is native to the southern part of the United 
States and to Mexico. It spread quickly in the eastern part of the United States 
and Canada in the 1940s after the release of captive specimens in the New York 
City area. The model is of the form 



N n+ i(x) = s(l- PA )N n (x) + (l-pj)F(N n (x))+ K A (\x-y\)s PA N n (y)dy 

J — OO 

(1.1) +/ Kj(\x-y\)pjF(N n (y))dy, n=l,2,— , 



where N n = J n + A n is the sum of the juvenile and adult bird densities in year 
n. It is assumed that juvenile birds born in one season mature by the next season. 
Adults survive from one season to the next with probability s. The number of 
surviving offspring born in year n is denoted by the fecundity function F(N n ). 
Juveniles and adults may differ in their probability to disperse as well as in their 
dispersal behavior if they disperse. The fractions of dispersing juveniles and adults 
are denoted by pj and p A , respectively. The dispersal behavior is given by the 
probability density functions Kj,K A , also called dispersal kernels. 

In this paper, we will always assume that the fecundity function is a monotone 
increasing and bounded function, for example a function of Beverton-Holt type or, 
equivalently, of Holling type II. The Beverton-Holt function originated as a stock- 
recruitment function in fisheries and is now very common in classic discrete-time 
population models. By appropriate rescaling, we can always assume that F is of 



Date: November 12, 2010. 

2010 Mathematics Subject Classification. Primary: 92D25; Secondary: 37N25; 39A22. 
Key words and phrases. Traveling waves; spreading speed; partially sedentary population; 
delay effect. 



1 



2 



THUC MANH LE, FRITHJOF LUTSCHER, AND NGUYEN VAN MINH 



the form 



(1.2) 



F(u): 



krMu 



M+{r-l)u 



with k + s = 1. All parameters are assumed positive. At low population density, the 
per capita number of offspring is fcr, and the carrying capacity of the non-spatial 
model N n+1 = sN n + F(N n ) is M. 

As noted in [25] . although this model was used to simulate the migration of 
house finches in [24], no mathematical analysis was given yet. A first attempt 
to study models with partially sedentary populations was made in |25) . These 
authors studied the special case that dispersal probability and dispersal behavior 
are independent of age-structure. More precisely, under the assumptions p,j = pa 
and Kj = Ka, the above model falls within the framework considered in [25]. It 
is the purpose of this paper to give a thorough study of the asymptotic behavior 
of solutions of Eq. (jl.ll) . (II .2[) with pa ^ p.j and Kj ^ Ka- Qur goal is to extend 
the fundamental theory for spreading speeds and traveling waves developed by 
Weinberger to this case. For the existence of traveling waves the fundamental 
assumption on compactness of [26l Theorem 6.6] is not satisfied for Eq. (jl.ip . nor 
the weak compactness condition in [TO]. We will use a different approach to prove 
the existence of traveling waves in this case. This idea is also employed in [12] . Our 
main result is Theorem 13 . 5 1 that complements some results in [25) . 

To simplify notation, we will assume that Kj = Ka = ■ K from here on. We 
return to the general case in the discussion. 

Notations and Assumptions. We denote by R the real line. We also denote 
by BAf(R,l) (_BC(R, R), respectively) the space of all measurable and essentially 
bounded real valued functions on R (the space of all bounded continuous real valued 
functions on R, respectively) with essential sup-norm. For a constant a we will 
denote the constant function R 3 x i— ► a by this number a for convenience if this 
does not cause any confusion. Cm stands for the set {/ S BC(R, M.)\f(x) S [0, M]}, 
and BM(M, [0, M]) := {/ e BM(R,R)\f{x) e [0, M]}. The metric on C M is defined 
by the sup norm. In _BA/(R,R) we use the natural order defined as u < v if and 
only if u{x) < v(x) for all x £ R. 

Unless otherwise stated, we assume that the parameters in the function F satisfy 
r > 1,M > 0. We also assume that K(\x\) is a probability density function defined 
on R and satisfies 



Most of the results in this section are derived from the general theory on spread- 
ing speeds in [30J . For later use in the paper we will discuss details of these results 
below. Let us define a dynamical system u n +i — Q[u n ] by setting 



Q[u]{x) = s(l - pa)u{x) + (1 - pj)F(u(x)) + K(\x - y\)sp A u{y)dy 



(1.3) 




for all n e R. 



2. Spreading Speed 




— oo 



(2.1) 
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for each u £ M(K,1). 

Lemma 2.1. Under the above notations and assumptions, the operator Q is an 
operator acting in £>M(R, R) leaving _BC(R, R) invariant with the following prop- 
erties 

i) Q[0] = 0, Q[M] = M, Q[a] > a /or all < a < M; 

ii) Ifu,ve BM (R, R) suc/i tfiai u > u, iften Q[u] > Q[«]; 

iii) Ifu n E £?C(R, R) smc/i £/ia£ u„ is convergent to u uniformly on each bounded 
subset o/M, i/ien <3[u n ](a;) is convergent to Q[u](x) for each x G R; 

iv) If a > M, then Q[a] < a; 

v) There is a constant 7 such that 7 < Qpy] < 7 /or a// 7 S [0,7), and 

Q[i\ = 7. 

Proof. Before proving the properties of Q we notice that from the definition of Q 
it maps BM(l,t) into itself. Next, to show that it leaves £>C(R,R) invariant it is 
sufficient to prove that the integrals 



(2.2) / K(\x-y\)u(y)dy = K{\t\)u(x 

J — OO J — OO 

OO /• 00 

(2.3) / tffla: - y|)F(u(y))dy = -/ #(|£|)F(u(z - 0)d£ 



depend continuously on x. To this purpose, for each e > we can find a large 
enough T > such that 

i-—T 



/— 1 poo 
K(\Z\)d£+ K(\S\)d£< 
-00 Jt 



4(H + ||F(h)||) 

On the interval [— T — 1,T + 1] the functions u and F(u) are both uniformly con- 
tinuous. Therefore, for a given xq there exists 1 > 8 > such that if \xo — x\ < 5, 
then, for all £ G [-T, T] 

(2.5) |u(s - £)-«(*«) - 01 < J 

(2.6) IFKxo-O)-^^-^))! <J. 
This yields that for all \xq — x\ < 8, 

/oo 
A-dClMas-Odel 
-oo 

■OO rOO 

<\l K(\m+ m\)d£\-\\u\\ 

T JT 

T K(mu(x ~o-u(x-om 

-T 

e e 

^ 7 + 7 < e - 
4 4 



Similarly 



That is the continuity of the integrals in (|2.2j) . proving that Q leaves BC(R,R) 
invariant. 



K(\Z\)F(u(x - 0)dC " / #(l£l)*Hs - OKI < e - 



4 



THUC MANH LE, FRITHJOF LUTSCHER, AND NGUYEN VAN MINH 



Now we prove the properties of Q: 

(i) : This property is clear because of the assumption s + k = 1. Indeed, Q[0] = 0, 
and 

krMM f°° 

Q[M]{x) = 3 (i- PA )M + (l-pj) M+ + J K{\x-y\)s PA Mdy 

f°° „n n krMM 
+ J_ oo K(\x-y\)pj M+ir _ 1)M dy 

= s(l - pa)M + (1 - pj)kM + sp A M + pjkM 
= (s + k)M 

(2.7) = M. 

Notice that F(a) > ka for all < a < M. Therefore, for all < a < M 

/oo 
K(\x - y\)spAOtdy 
-oo 

K(\x — y\)p,jkady 

= s(l — pa)ch + (1 — Pj)ka + sp^a + pjka 
— (s + k)a 

(2.8) = a. 

(ii) : This property follows since the function F is increasing. 

(iii) : This property is clear from the definition of Q. 

(iv) : We have 

. r . , kr(M + (r - l)u) - kru(r - 1) , 

b (it) — k = 7— — -, , ,„ k 

w (M + (r - l)u) 2 

krM 



(M + (r - 1» 2 

Notice that F'(u) - fe < if u > M. Therefore, F(a) < ka if a > M. This yields 
that if a > M 

/oo 
K(\x — y\)spAOtdy 
-oo 

K{\x - y\)pjkady 

(2.9) = a. 

(v): Actually we can choose 7 = M. In fact, for each 7 g (0, M) we have ^(7) > 
k-f. Therefore, by the above computation (12.81) . 7 < (3(7) < Q[M] — M for all 
j€[0,M}. ' ' □ 

Next, we apply the general theory on spreading speeds in 26 to our dynamical 
system 

(2.10) u n+ i = Q[u n ], n=l,2,---. 

Basically, according to the theory of for our model the following procedure of 
defining the spreading speed is valid: let us choose a function ip : K — > M such that 

i) if is continuous and non-increasing, 

ii) ip(-oo) := limt^-oo <p(t) £ (0,M), 
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iii) ip(s) — for s > 0. 
We then define an operator R c [-] on the space Cm for every constant c as 

(2.11) R c [u](s) := max{ip(s),Q[u(c + ■)](&)}, set, 
and a sequence of functions {a„(c; •)} by 

(2.12) a n+ i := R c [a n ], a = tp. 

As shown in [26] the sequence {o n (c; •)} is increasing and bounded, so for each 
s £ R, we obtain the pointwise limit 

(2.13) lim a„(c; s) = a(c; s) , set. 

n— >oo 

Obviously, < a(c; s) < M for all sgl. The following number is called spreading 
speed for our model 

(2.14) c* := sup{e|a(c; +oo) = M}. 

Applying the theory in [26, Lemma 5.2 and Proposition 5.1] to our model gives the 
following: 

Lemma 2.2. i) For each eel 

(2.15) a(c;-oo)=M; 
ii) If c > c* , t/ie« 

(2.16) a(c; +oo) = 0. 

Weinberger [26] proved that if there exists a bounded non-negative measure 
777,(2;, dx) on K such that 



(2.17) Q[u](x) < / u(x- y)m(y, dy), u£C M , 

J — 00 

then 

1 r 00 

(2.18) c* < inf - In / e^mtx, dx). 

And if there exists a bounded non- negative measure l{x, dx) with property that 
J_ l(x, dx) > 1 and Q[u](x) > J u{x — y)l{y, dy) for all u such that < u(x) < e, 
then 



(2.19) 



1 r 00 

c* > inf - In / e^ltx, dx). 



Below we follow the argument of [33] to give an estimate of the spreading speed. 
We see that since for u £ Cm F{u) < kru, 
(2.20) 

/OO 
K (\ x - V\) (sp A +pjkr) u{y)dy. 
-00 

If we let 

m(a;, dx) = [s(l - pa) + (1 - pj)fcr] 80 + K(\x\) (spa + pjkr) , 
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where So is the Dirac delta measure, then (|2.17|) holds. On the other hand, for each 
1 < n < r there exists e > such that F(u) > kr\u for u such that < u(x) < e 
for all Therefore, for such u 

(2.21) 

p OO 

Q[u}(x)>[s(l-p A ) + (l-p J )r 1 ]u(x) + / K(\x-y\)(sp A +p J r 1 )u(y)dy. 



We let 

l(x, dx) := [s(l - p A ) + (1 - 2v)feri] <5o + K(\x\)(sp A + pjkn). 

Then, X!^i(x, doc) > 1 and (l2~T9l) holds. Next, since r% can be chosen arbitrarily 
between 1 and r, we have 

inf - r 



1 f 00 

(2.22) c* = inf - / (sp A +pjkr)e^ x K(\x\)dx+[s(l~p A ) + (l-pj)kr}. 



3. Traveling Waves 

This section contains our main results of the paper. 

Definition 3.1. A monotone traveling wave solution with speed c connecting to 
M, or for short a traveling wave, of Eq. is defined to be a non-increasing 

continuous function w such that lim^^oo w(x) = 0, lim 2 ,_ > ._ 00 w(x) — M, and 
N n (x) := to (x — nc) is a solution of Eq. (jl.lj) . 

If we substitute N n into (|l.ll) . we will have 
u>(x — (n + l)c) = s(l — p_4)io(a; — nc) + (1 — pj)F(w(x — nc)) 

K(\x - y\)sp A w(y - nc)dy 



(3.1) +/ K(\x-y\)pjF(w(y-nc))dy, n=l,2,---. 

If we set £ := a; — (n + l)c and z := y — nc then, the above equation becomes what 
is called "wave equation" associated with (|1.1[) : 



w(0 = s(l-p,iM£ + e) + (l-jv)F(u;(£ + c))+ / + c - z\)sp A w{z)dz 

J —OO 

rOO 

(3.2) +/ K(\£ + c-z\) Pj F(w(z))dz, n = l,2,---. 



For each u G Cj\/, or more generally BM(K, K), and celwe set 
(3.3) B c [u]{x) = s(l-p A )u(x + c) + (l-pj)F(u(x + c)) 

/•OO 

C c [u)(x) = / + c - y\)sp A u(y)dy 



(3.4) +/ A-(|i + c-i/|)pjF(u(i/))di/. 

^ — oo 

Obviously, a monotone traveling wave solution to (jl.ll) is a non-increasing continu- 
ous function w with w(— oo) = M, w(oo) — that is a fixed point of Q c := B c + C c . 
We remark that although there are several extensions of Weinberger's theory on 
the existence of monotone traveling waves, to our best knowledge the existence of 
monotone traveling wave to (|1.1|) is still open. The reason is the operator Q c does 
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not satisfy any compactness conditions listed in [26] , [10] . Volkov and Lui extended 
Weinberger's theory to a class of systems without compactness conditions. How- 
ever, the model considered in [2 5) includes our model only if pa = pj, that 
means that age structure does not affect migration behavior. However, it is more 
realistic to make the assumption that pj ^ pa ■ This assumption makes the problem 
of studying the existence of monotone traveling waves much harder. Below we will 
prove the existence of monotone traveling waves of under assumption that the 
function F in (fTTT]) is of the form ([L2")l . 

Lemma 3.2. Assume that 

(3.5) s{l- PA ) + {l-pj)kr < 1. 

Then, for each given w G Cm (w G BM(R, [0, M]), respectively) the equation oper- 
ator u — B c [u] = w has an unique solution u in Cm (in BM (R, [0, M\), respectively) 
which will be denoted by u :— G c w. 

Proof. Consider the function 

(3.6) g(x) = 8(1 - p A )x + (1 - PJ ) * TMX xg[0,M]. 

M + (r — l)x 

For all x G [0,M], 

, . , . . . krM(M + (r - l)x) - krMx(r - I) 
g\x) = s(l - p A ) + (1 - pj) ^ r \, f | /_ ' 2 K '- 



= s(l - p A ) + (I - pj) 
< s(l-p A ) + (l-pj) 



(M + (r - l)xf 
krM 2 



(M + (r- l)x) 2 
krM 2 



(M + (r - 1) ■ 0) 2 
= s(l - P a) + (1 -pj)kr 
< 1. 

Therefore, there exists a positive < p < 1 such that 

(3.7) 0< sup \g'(x)\<p. 

XE[0,M] 

Next, we solve the equation u — B c [u] = w for each given w G C[0, M]. Note that 
in this case, B c is a strict contraction because 

||-B c [wi] - B c [u 2 }\\ = sup \g(ui(x + c)) - g(u 2 (x + c))| 
< sup \g'(i)\ ■ \u\(x + c) - U2(x + c)\ 

(3.8) = p[|«i-u a ||. 

Therefore, by a standard argument we can prove the existence of (/ — Bc)^ 1 that 
is Lipschitz continuous. □ 

Lemma 3.3. Assume that the kernel K(\x\) satisfies the above mentioned condi- 
tions. Then, for each monotonous u the function C c [u] is continuous. 

Proof. The proof can be done in the same way as in that of [HI Lemma 3.5]. □ 



s 
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Remark 3.4. We notice that although the operator B c is a strict contraction in the 
uniform convergent topology it is not a strict contraction in the norm of compact 
open topology as defined in [TU]. Moreover, operator C — c is not compact the 
norm of compact open topology because the kernel K may not be continuous, so 
condition (A3) in [10] is not satisfied with the operator B c + C c . That is, the theory 
of traveling waves in [2 6) as well as its extension in [lOj does not apply to this case. 

The following is the main result of the paper: 

Theorem 3.5. Let all assumptions in Lemmas \3.S\ and \3.3\ be satisfied. Then, if 
c > c* , then there exists a monotone traveling wave to Eq. M.l)) . 

Proof. Let the function a(c; ■) be defined as in (12.131) . Set (f)%(s) :— a(c;s). Note 
that since a(c; •) is non-increasing and bounded it is a measurable and bounded 
function on R. We define a sequence 

(3.9) 4>n+i = Q C [<M, n>l, n=l,2,---. 

We now show that {<fi n } is a non-increasing sequence in BM(K, [0, M]). In fact, by 
definition of the sequence {a„(c; •)} we have 

a n+1 (c;s) := max{^(s), Q[a n (c; ■ + s + c)](0)} 

(3.10) = max{</J i (s),Q[a„(c;- + c)](s)} 

> Q[a n (c;- + c)](s), 

(3.11) = Q c [a n {c;-)}(s) 
Therefore, 

(3.12) a(c;s) > Q c [a(c; .)](*). 
That is, 

(3.13) 4>x > 02. 

Since Q c is order-preserving, by the definition of 4> n +i, (|3.13[) yields that 

(3.14) <f> n > <f> n+l , n e N . 

Therefore, the sequence {4> n } is pointwise non-increasing and bounded below by 
zero (because these functions are non- negative) , so it has a limit W that is a non- 
increasing function, so it is measurable. We will show that W is a traveling wave 
solution to Eq. ([TTT]) . By Lemma l2~2l W(-oo) = M. Next, < ^(+00) < 
(f>i (+00) = 0, so, W(+oo) = 0. In particular, W is a fixed point of Q Cl that is 

(3.15) W = Q C [W] := S C [W1 + C C [W]. 

We need only to show that W is continuous for it to be a traveling wave as in our 
definition. Since Q3.15P is equivalent to the following 

(3.16) W - B C [W] = C C [W], 
by Lemma I3T21 it is equivalent to 

(3.17) W=[I-B C ]- 1 C C [W], 

By Lemma l3~3l the function C C [W] is continuous. In turn, by Lemma l3?2| G C [W] = 
[I — B C ]~ 1 C C [W] is continuous. Therefore, W is continuous, and thus it is a wave 
solution of Eq. (fTTTj) . □ 
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Remark 3.6. As in 26, TQl [12] by using Lemma l2~2l we can easily show that if c < c* 
the traveling waves do not exist. With this said the spreading speed c* is exactly 
the minimal wave speed of traveling wave solution to Eq. (11.11) . 



In Section 3 we considered the existence of traveling waves when F is of the 
form (|1.2[) . The function F can chosen to be a more general one that satisfies the 
following conditions with given positive numbers r, fc, M: 

(HI) F e C l [0,M\; 

(H2) F(0) = 0, and F(M) = kM, s + k = 1; 

(H3) F(u) > ku, for u e (0,M); 

(H4) F'(u) > 0, and .F'(O) = fcr, where r > 1 is a given constant; 

(H5) F(u) < kru, for u 6 [0,M]; 

(H6) F'(x) is non-increasing on [0, M]. 

Then, the statement of the main results as well as its proof are unchaged. 

Results of the previous sections can be easily extended to the case the habitat 
is multiple dimensional K d with d = 2,3,.... There are no big changes in the 
statements of the results. And the ideas of proofs remain similar. 

We turn to the more biological aspects of our work. Condition (|3.5|) imposes 
an upper bound on the per capita offspring production or, more precisely, on the 
proportion of individuals who do not disperse. The original model by Veit and Lewis 
[24] allowed for different dispersal behavior of juveniles and adults. The theory 
presented here easily extends to the case K A ^ Kj. Most importantly, Lemma [3~3l 
holds if the continuity conditions holds for both kernels. Formula (|2.22|) for the 
spreading speed becomes c* = inf M> o — /t(/i), where 



(4.1) k(/i) = / [sp A K A (\x\) +pjkrKj(\x\)]e^ x dx + [s(l - p A ) + (1 - pj)kr 



More important differences between the model by Veit and Lewis and our analysis 
here is that they considered the function F to describe a strong Allee effect and 
the dispersal probabilities, pj,p A , to depend on population density. A strong Alice 
effect occurs if the per capita population growth rate is highest for intermediate 
population densities so that the population actually declines for small densities. 
While the existence of a spreading speed in the presence of an Allee effect is still 
guaranteed by Weinberger's theory, the existence of traveling waves for that case 
is a wide open question. Equally open is the question of traveling waves and even 
the existence of a spreading speed for models with density-dependent dispersal 
probability. Some preliminary results and caveats were obtained in [llj . These 
questions remain the subject of our future investigation. 
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